Radiation content of Conformally flat initial data 
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We study the radiation of energy and linear momentum emitted to infinity by the headon collision 
of binary black holes, starting from rest at a finite initial separation, in the extreme mass ratio limit. 
For these configurations we identify the radiation produced by the initially conformally flat choice of 
the three geometry. This identification suggests that the radiated energy and momentum of headon 
collisions will not be dominated by the details of the initial data for evolution of holes from initial 
proper separations Lo > 7M. For non-headon orbits, where the amount of radiation is orders of 
magnitude larger, the conformally flat initial data may provide a relative even better approximation. 
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I. INTRODUCTION 

As we get closer to the first direct detection of grav- 
itational waves, there is a renewed interest in the con- 
struction of astrophysically relevant initial data for bi- 
nary black holes, a priori thought to be the most violent 
source of gravitational radiation in the Universe. The 
traditional choice for initial data had been the classical 
work of Bowen and York 0. More recently, attention 
has been directed toward 'thin-sandwich' 2, 3] and Kerr- 
Schild ansdtze 0, 0| initial data, and comparisons with 
the Bowen- York choice have been carried out on the ini- 
tial shces IE S ■ Post-Newtonian hmit inspired data 
are also beginning to appear in the literature 0, . 

The most important question about the choice of ini- 
tial data is the effect of that choice on the gravita- 
tional waveform at infinity. Answering that question is 
very difficult if the binary black holes have comparable 
mass [III, H^i HJi HH I but is relatively straightforward 
in the particle limit. For this reason, we have previously 
studied the headon collision of binary black holes in the 
extreme mass ratio regime|0, ^3 

hj. This allowed us 
to compare 'on slice' initial data and truly evolved data. 
In It', we focused particular attention on the im- 
portance of the choice of the extrinsic curvature for the 
initial data, and we considered only the usual choice for 
the initial three metric: conformal flatness. In this re- 
port we turn to the question omitted from |0, , and 
we look briefly at the importance of conformal flatness. 
We consider time symmetric initial data so that the nat- 
ural choice for the extrinsic curvature is for it to vanish. 
We use the Brill-Lindquist type of conformally flat initial 
data. (We have shown in Ref. J/^] that the total radiated 
energy is not very sensitive to this choice.) We evolve the 
time symmetric initial data, compute the radiated wave- 
forms, and try to identify what part of the radiation can 
be ascribed to the conformally flat choice of the initial 
three metric. We do this differently for large and for 
small initial separation of the holes. In the case of large 
initial separation, we can identify an early feature of the 
waveform that is clearly produced by the initial data. For 



the case of small initial separation, we introduce a more 
speculative measure of the radiation ascribed to the con- 
formally flat initial data; we consider it to be the excess 
radiation above the minimum of a one-parameter family 
of initial data choices. 



II. RESULTS 

We consider a particle of mass mo, to be a first-order 
perturbation on the background of a Schwarzschild hole 
of mass M . The particle is initially at rest at a proper 
distance Lq from the horizon of the Schwarzschild hole. 
For the radial infall of the particle, the odd parity per- 
turbations of the spacetime vanish. For the even parity 
perturbations we use the formulae and numerical tech- 
niques described in ^3 to integrate the Moncrief-Zerilli 
equation 

[dl. ~ d1 - V,(r)\ v., = S,{T^{t), r) , (2.1) 

for the even parity wave function ipi{t,r) for each £-pole 
mode of the field. 



A. radiated energies 

The radiated even-parity energy at infinity is computed 
from 

-=l^n^V^^(a.^.(M))^ (2.2) 

The computed energy Eg for i—2. . . 5, reported in TablelH 
and Fig. ^ show the expected exponential decrease of 
energy with increasing £0, 0. The plot also displays 
an unexpected feature: For small separations there is a 
rise in the radiated energy. The relative importance of 
this feature grows with increasing £, and for £ > 3 the 
feature is so strong that radiation from small Lq/AI is 
greater than that for infall from infinity. 
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FIG. 1: The contribution to the radiated energy at infinity of 
the first four multipoles, as a function of the initial location 
of the particle falling from rest. 



This increase in Ei for small Lg /M is certainly related 
to the choice of initial data. As a first step in under- 
standing this relationship, in Fig. |21we display dip2/dt, 
the quantity that dominates the energy in Table ^ and 
Fig. n The waveform of d4>2/dt, for an observer at 
r* /2M = 100, shows the characteristic excitation and 
ringdown of quasinormal oscillations. But at early times, 
t/2M ~ 100, a small feature appears representing the 
evolution of the initial disturbance of the Schwarzschild 
background around the particle position. We have iso- 
lated this small feature and have computed its contri- 
bution SE2 to the quadrupole radiated energy E2. The 
values of dE2 are listed in TableHl and the circles in Fig.|3| 
show the fractional energy 6E2/E2 as a function of the 
initial separation. The circles in Fig. |31 show a marked 
increase in the fractional energy as Lq/M decreases, but 
this method of ascribing energy to the initial data can 
only be extended down to Lq/M ~ 10. For Lq/M < 10 
the small early feature in the waveform cannot be cleanly 
distinguished from the initial excitation of ringing. (The 
error in dE2 reaches several percent for the values tabu- 
lated in Table HI) 

One can make a speculative estimate of the small- 
Lq initial-value energy by extrapolation. A best fit of 
the form SE2/E2 — a(Lo/Af)^ requires a = 7.26 and 
b = -2.135. This fitted curve reaches SE2/E2 ~ 0.5 for 
Lo/M « 3.5. 

There is a completely different way to estimate the en- 
ergy associated with the conformally flat initial data. To 
do this, we follow the approach of Martel and Poisson|23 . 
In the Regge- Wheeler |2lj gauge for the multipole pertur- 
bations of Schwarzschild spacetime, the only nonvanish- 
ing even parity perturbations on the t = constant initial 
hypersurface are ijf'" and X^™. The initial value equa- 
tions contain a free functional degree of freedom in 



FIG. 2: Waveforms of dtp2/dt at r* /2M = 100. Note the 
small feature located around t/2M = 100, associated with 
the initial data radiation content. 
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FIG. 3: The estimated fraction of radiated £ = 2 energy that 
is due to the choice of conformally flat initial data. The cir- 
cles, for Lo/M > 9.56 represent fractional energy in the early 
feature in the waveform. The squares, for Lo/M < 8, are 
based on the energy in excess of a baseline model of initial 
data that is not conformally flat. 



and K^"^. The usual choice, i/l™ — K^™, corresponds to 
a conformally fiat perturbed three metric. Here, instead, 
we consider the one parameter family 



= aK' 



(2.3) 



In Rcf. 20J, the energy radiated has been studied as a 
function of a, and the value of a has been found for 
which the radiated energy is minimum. This turns out 
not to be the conformally flat choice a = 1. (See Fig. 12 
in Ref. ^20].) For any value of Lq/M, we take the a- 
model of Eq. (|2.3|) with the minimum radiated energy as 



TABLE I; Energy radiated in units of ml/2M 



Lq/M E2 



E. 



E. 



Total 



104.3 

63.8 

33.0 

22.6 

17.3 

14.0 

12.9 

11.8 

10.7 

9.56 

7.19 

5.94 

5.28 

4.59 

4.23 

3.86 

3.05 

2.30 

1.29 



0.0179 

0.0174 

0.0164 

0.0157 

0.0151 

0.0147 

0.0145 

0.0143 

0.0142 

0.0141 

0.0140 

0.0154 

0.0157 

0.0150 

0.0140 

0.0125 

0.00824 

0.00424 

0.00095 



0.00214 
0.00207 
0.00194 
0.00185 
0.00176 
0.00169 
0.00166 
0.00165 
0.00162 
0.00160 
0.00170 
0.00164 
0.00181 
0.00223 
0.00244 
0.00256 
0.00219 
0.00122 
0.00020 



0.000318 
0.000289 
0.000262 
0.000257 
0.000245 
0.000234 
0.000237 
0.000234 
0.000232 
0.000237 
0.000271 
0.000336 
0.000350 
0.000378 
0.000429 
0.000505 
0.000592 
0.000383 
5.24E-05 



9.69E-05 
5.61E-05 
4.20E-05 
5.31E-05 
5.30E-05 
4.92E-05 
5.28E-05 
5.20E-05 
5.39E-05 
5.99E-05 
8.11E-05 
0.000105 
0.000135 
0.000145 
0.000144 
0.000150 
0.000189 
0.000140 
1.64E-05 



0.0205 

0.0198 

0.0187 

0.0179 

0.0172 

0.0167 

0.0165 

0.0162 

0.0161 

0.0160 

0.0160 

0.0174 

0.0180 

0.0177 

0.0170 

0.0157 

0.0112 

0.00595 

0.00121 



.49226- 

.36956- 

.30626- 

.77076- 

.15366- 

.25576- 

.33446- 

.38116- 

.45496- 

.57306- 

0.140 

0.362 

0.465 

0.509 

0.494 

0.485 

0.409 

0.354 



r 



the "baseline" model. We subtract the radiated energy 
for this baseline model from the radiated energy for a 
conformally flat model and consider the excess energy to 
be an artifact of the conformally flat initial data, radia- 
tion energy that is, in a sense, contained within the initial 
data. This excess energy, for ^ = 2 and Lq/M < 8, is in- 
cluded in Fig. 121 It may be interesting that the maximum 
of this excess energy is located near a proper separation 
Lq/M « 4.5, in rough agreement with the previous ex- 
trapolation. 



B. radiated linear momentum 

An astrophysically interesting quantity to compute is 
the radiated linear momentum and the corresponding re- 
coil velocity of the final black hole. If the recoil velocity 
can be comparable to the escape velocity of a typical 
galaxy (100 - 1000 km/s), or a cluster (10 - 20km/s), 
recoil-driven escape can affect the growth rate of mas- 
sive black holes in the core of such star systems. The 
computation of this quantity for a particle released from 
rest at infinity in a Kerr background was performed in 
Refs. |23,|2^. For collisions of comparable mass holes the 
radiated linear momentum in a headon collision was com- 
puted in Ref. 24] in the close limit, and for several larger 
initial separations in full 2D nonlinear general relativity 
in Ref. [23 • We present here explicit results of the ex- 
treme mass ratio limit in order to provide benchmarking 
data for future full nonlinear numerical work. 

The radiated linear momentum along the axis of the 




FIG. 4: Recoil velocity of the system for different initial 
proper separations of the holes. The dashed line indicates 
an estimate of recoil velocity for a baseline model of initial 
data that is not conformally flat. 



collision is computed as the correlation of two successive 
multipole contributions |2^ 



dt 



= lim 



1 

32^^ 

i=2 



{£ + Sy.dtMr,t)dt^e+i{r,t) 
{£-2)1 ^{2£+l){2£ + 3) 



(2.4) 

Figure0]displays the results of computing the recoil ve- 
locity of the binary system by summing over i ~ 2 ■ ■ ■ 5. 
A notable feature here is that radiated momentum is a 
maximum for a particle that starts at rest from Lq /M ^ 
4.5, not from infinity. For this momentum maximum, 
the recoil velocity reaches almost 400(r7io/Af)^ km/s com- 
pared to 250 (mo /Af)^ km/s for a particle from infinity 
(See Table As we did with the radiated energy, we 
can make a speculative estimate of the fraction of the re- 
coil velocity that may be ascribed to the conformally flat 
initial data. We take 5v to be the excess recoil velocity 
above that of the baseline (energy minimum) model. At 
the peak recoil velocity we find that 5v is 1.0 times the 
recoil velocity for the baseline model. Thus the part of 
the recoil that we view as an artifact of the initial data 
is as big as the baseline value. This effect is more pro- 
nounced than that for energy because the radiation of 
momentum involves ^ = 3 as much as it does £ — 2; as 
shown in Fig.^ the sensitivity of dipg/dt to initial data 
increases with increasing £. 



III. DISCUSSION 

To assess the importance of the choice of the initial 
three geometry on the computation of radiation from 
black hole mergers, we have used the very simplest model: 
the extreme mass limit, a nonrotating hole, and radial 
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TABLE II: Linear momentum radiated 





Lo/M 


( M lmn\^it km/s 


100 


104.3 


249 


60 


63.8 


248 


30 


33.0 


237 


20 


22.6 


230 


15 


17.3 


222 


12 


14.0 


218 


11 


12.9 


217 


10 


11.8 


216 


9 


10.7 


219 


8 


9.56 


221 


6 


7.19 


257 


5 


5.94 


320 


4.5 


5.28 


352 


4 


4.59 


372 


3.75 


4.23 


371 


3.5 


3.86 


358 


3 


3.05 


271 


2.6 


2.30 


144 


2.2 


1.29 


24 



infall from rest. Our results suggest that for proper dis- 
tance Lo larger than around 7M the energy and momen- 
tum radiated are probably not significantly contaminated 
by the initial choice. In the astrophysically more interest- 
ing case of inspiralling binary black holes in quasicircu- 
lar orbits, we expect the radiation will be affected much 
more by the choice of the extrinsic curvature, than the 
choice of the three metric. In that case, the conformally 
flat choice of the three-geometry seems adequate for most 
of the applications. A study of the initial data choices 
for quasicircular orbits in the small mass ratio limit is 
currently underway by the authors. 
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